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ABSTRACT 

In the inner regions of an accretion disk around a black hole, the presence 
of non-thermal protons would via p-p collisions produce electrons, positrons and 
^ . 7-rays. For such a scenario, the steady state electron-positron distribution is 
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computed taking into account Compton cooling, e~e+ pair production (due to 
7 — 7 interactions) and pair annihilation. 

The resultant spectra has a broad feature around 1 — 10 MeV which may be 



on the non-thermal proton fraction can be obtained, which may put constraints 
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00 ! tested with observations by INTEGRAL and GLAST. Using the available OSSE 
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data for the black hole system, GRS 1915+115 we illustrate that an upper limit 

O . 

. on the acceleration process active in this system. 
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d • 1. Introduction 

^> : 

^ I Black hole X-ray binary systems are observed primarily in one of the two long term 

^ I spectral states : the hard-state and the soft-state. In the so called hard state, the spectrum 
of the system can be described approximately as a hard power-law with a spectral index 
r 1.7 and a cutoff around 100 keV, which suggest a thermal Comptonization origin 
(Shapiro et al. 1976; Liang & Price 1977). Indeed, detailed fits to broad band X-ray spectra 
confirm that the spectra is well represented by Comptonization of soft photons by a plasma 
having temperature T 50 keV and optical depths of order unity (Gierlinski et al. 1997). 
In contrast, the soft spectral state consists of a black-body like component (kT 1 keV) 
which generally dominates the luminosity, and a power-law tail with a photon index F ^ 2.5 
with no detectable cutoff up to ^ 600 keV (Zdziarski et al. 2001). The steepness of the 
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spectra and the absence of a spectral break indicates that this component arises due to 
Comptonization of soft photons by a non-thermal electron distribution. Indeed, the spectra 
can be well fit by a model where soft photons are Comptonized by a plasma consisting of both 
thermal and non-thermal electron distributions in a "hybrid" state (Zdziarski et al. 2001; 
Poutanen & Coppi 1998). Comptonization by thermal electrons {kT ^ 1 keV) gives rise to 
the soft component, while the power- law tail is produced by the non-thermal electrons. In 
this model, the steady state non-thermal electron distribution is computed self-consistently 
by assuming that there is an injection of non-thermal particles into the system where they 
are cooled by Comptonization and Coulomb interactions (Gierlinski et al. 1999). The salient 
results of this analysis suggests that the thermal and non-thermal components of the electron 
distribution co-cxist in the same region and that the power in the non-thermal injection rates 
may be ~ 10 — 20% of the luminosity of the source. 

The hard X-ray emitting region in the soft state can possibly be a corona or a flare on 
top of the cold accretion disk (Haardt & Maraschi 1993; Poutanen & Fabian 1999). While 
the energy released in a corona/flare is probably related to magnetic field reconnections, the 
details of the process are largely unknown. It is also not clear what kind of an electron accel- 
eration mechanism (which is inferred from the non-thermal spectra) is active in these sources. 
If this mechanism is mass independent then one expects protons also to be accelerated, to 
form a non-thermal proton distribution. Protons may also be independently accelerated by 
scattering off magnetic "kinks" in a Keplerian disk (Subramanian et al. 1996). Thus ev- 
idence for a non-thermal proton distribution is expected to give valuable insight into the 
nature of the accelerating processes, and hence constrain the nature of the X-ray emitting 
region. Moreover, these protons may also escape from the inner regions and form relativistic 
jets (Subramanian et al. 1999) observed in some of these sources. 

Relativistic protons can undergo inelastic collisions with thermal protons to produce 
charged and neutral pions which, in turn, decay into electrons, positrons and 7-rays. Since 
these electrons and positrons are primarily produced at high energies (7e ~ 200), their steady 
state distributions are expected to be steep (typically oc 7^). Thus their photon spectral sig- 
nature will be distinct from the non-thermal electron distribution produced directly by the 
acceleration process (typically oc 7~^). Pair cascades initiated by injection of pairs in the 
system e.g. (Lightman & Zdziarski 1987; Svensson 1987) and pair production due p — p 
reactions (Stern et al. 1992) have been studied extensively. In these works, Monte Carlo 
simulations have been undertaken to compute correctly the emergent spectra from a pair 
dominated plasma. Our motivation is to compute these spectral signatures for parameters 
pertaining to black hole binaries in the soft state which can be used to constrain the fraction 
of non-thermal protons in the system. We limit the study to plasmas which are not pair 
dominated which allows for approximations that make the analysis relatively simpler. An 
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important effect, considered here, is the pair production of 7-ray photons in interaction with 
copious soft X-ray photons of energy ~ 3 keV that exist during the soft state. In this respect, 
this work differs from previous computation of spectra from non-thermal electron distribu- 
tions due to pion decay , where pair production due to the presence of copious soft photons 
are not considered (Eilek & Kafatos 1983; Mahadevan et al. 1997; Markoff et al. 1999) and 
those where the soft radiation were assumed to be UV photons (Stern et al. 1992), scenarios 
more relevant to under- luminous black hole systems and AGN. 

In the next section we describe the model and the assumptions made to compute the 
steady state non-thermal electron/positron distributions and the resultant photon spectra. 
In §3, generic results of the computation are presented along with a specific application to the 
black hole system GRS 1915-1-105. The work is summarized and main results are discussed 
in the last section. 



2. Steady state electron/positron distribution and spectra 

We consider a uniform sphere of non-relativistic thermal plasma with number density, 
TiT and radius R. To represent the soft spectral component, it is assumed that inside the 
plasma, there is a Wein peak photon density n^,s{e) at temperature Tg (where e = hv/rrieC^ 
is the normalized photon energy) and normalization determined such that the photon energy 
density 

/•oo 

Ury^s — T^ef? / en^^s{^)de — Ls{l -\- t) /{cAt^I^) (1) 
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Here Lg is the luminosity of the soft component and r = nrpcrTR is the Thompson optical 
depth. Since the acceleration mechanism is unknown, we assume that a small fraction / of 
the thermal protons form a non-thermal power-law distribution i.e 

cx — 1 'y 

nNxilp) = friT ( — ^)~" for 7p > 7p,^i„ (2) 

where 7p is the Lorentz factor of a proton and ^p^rnin is the minimum 7p of the distribution. 
It has been assumed here that the particle index a > 2 and that the maximum Lorentz 

factor, ^p^max -^'^ lfp,min- 

These non-thermal protons interact with the thermal protons via p — p reactions and 
subsequently produce electron/positron and 7-rays by pion decay. The threshold energy 
for such p — p reactions is ^p^thres ~ 1.5 while the other possible reactions ( e.g. p — 7 
pair production) have a much higher threshold '~fp,thres ~ 300 and are neglected. It is not 
clear to what maximum Lorentz factor ^p^max would the unknown mechanism accelerate the 
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protons. The absence of a spectral signature due to such reactions may simply indicate a 
lower value of jmax ^-nd not necessarily mean the absence of an efficient proton acceleration 
mechanism. Moreover, such spectral signatures would be sensitive to the proton power-law 
index, a assumed in this work, while the results obtained below are fairly independent of 
such details of the acceleration mechanism. Nevertheless for typical parameters used in this 
work, we a posteriori check the pair production rate due to the p — 7 reaction and confirm 
that it may be neglected. 

The steady state positron particle density, ne+{j), is obtained by solving the kinetic 
equation 

d 

^[(7c + 7comp)^e+(7)] + ne+{^) = Q+,'rM + <3+,pp(7) (3) 

while the electron particle density, n^,- (7) , is obtained from 

d 

^[(7c + 7comp)rie-(7)] = Q-jM + Q-,Pp{l) (4) 

Following (Eilek & Kafatos 1983), e+e" production rate due to the p — p process is given by 

Qe±,pp(7e) =nTc\ — 3/4^ 3/2.n/2 

where 

7p,/.(7e) = 1 + [7*7e + (7*' - - 1)^/' - 1]'/=^ (6) 

7.,/(7e) = 1 + [7.7e - (7.^ - - ^fl' - ^f' (7) 

where 7^ = 70. The approximate but analytical cross-section ( Ot^±) is tabulated for different 
energy ranges by (Eilek & Kafatos 1983) 

Pair production rate from photon-photon (i.e. 77) interaction is approximated to be 

Q±,77(7e) ^ y ^7(276 - e)n^(e)a^^(e, 27 - e)de (8) 

where it has been assumed that for two photons annihilating with energies e and e' the 
resultant e~ or e"*" Lorentz factor is f» (e-|-e')/2. The approximate form for the cross-section 
(7^^(e, e') is given by (Coppi & Blandford 1990). The positrons primarily annihilate with the 
background thermal electrons at a rate given by 

ne+(7e+) = ne+(7e+)nT(Te+e-(l,7e+)c (9) 



where the approximate form for the cross-section is used (Coppi & Blandford 1990). The 
annihilation of electrons has been neglected in eqn (4). 7c and 7comp ^-re the rates of change 
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of Lorentz factor due to Coulomb and inverse Compton cooling respectively. The inverse 
Compton cooling is primarily due to the up-scattering of soft (~ 3 keV) photons to higher 
energies by non-thermal particles. Since in general the leptonic Lorentz factor may be 
> 100, the scattering in the rest frame may not take place in the Thompson limit and the 
general expression for the inverse Compton photon production rate and cooling (Jones 1968; 
Blumenthal & Gould 1970) have been used. 

The equilibrium photon density inside the sphere is a solution of 

Q^,IC + Q7,PP + Ql,e+e- = n^(e)[i?77 + ---] (10) 

K[i + rKN[e)) 

where the rate of photon annihilation is given by 

R-y-y{e) — c J nj{e')aj-y{e,e')de' (11) 
with the cross-section approximated to be (Coppi & Blandford 1990), 
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1)'/% 1-1/2 , 3 



a^^{e, e') = ar^ ^7^(7;^" ^ + 7ln(a;))H(a; - 1) (12) 



X 
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Q-yjc, Q'y,pp Q'y,e+e- the photou productiou rates due to inverse Compton, p — p 
interaction and pair annihilation, respectively. The latter is assumed to be 

Qj,e+e-{e) = 4nTne+(2e - l)a^+^-{l,2e - l)c (13) 

The photon production due to p — p reaction is taken to be (Eilek & Kafatos 1983) 



m 



where 7^^ = 7^ — 1 and 



The approximate but analytical cross-section ( a-„o) is tabulated for different energy ranges 
by Kafatos and Eilek (1983). The last term in eqn.(lO) represents photons escaping from 
the system and tkn is the optical depth taking into account the Klien-Nishina cross-section. 

Although the Thompson scattering depth r of the thermal electrons could be of order 
unity, the corresponding Khen-Nishina optical depth for high energy (> 1 MeV) photons is 
small. Hence scattering of high energy photons with thermal electrons has been neglected. 
It should be emphasized that the above equations are valid only when the number density of 
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non-thermal leptons is much smaller than the thermal one ut- This assumption is violated 
when the plasma is pair dominated. It is known, that for some regions of parameter space, 
a pair runaway (or cascade) may be initiated and no steady state solutions for the leptonic 
density will exist. Solutions close in parameter space to this pair run away limit, may indeed 
be pair dominated and hence will not be accurately described by the above formalism. 
However, pair dominated plasma naturally exhibit a strong annihilation line at 511 keV 
which is not generally observed in black hole binaries. 

Equations (10), (3) and (4) are solved iteratively to obtain the steady state particle 
densities and radiative flux as a function of six parameters: the Thompson optical depth r, 
the size R, the soft photon luminosity Lg, the soft photon temperature T^, the fraction of 
non-thermal protons / and the non-thermal proton index a. 



3. Results 

The computed spectra for different values of the non-thermal proton fraction / and 
for typical galactic black hole parameters are shown in Figure 1 while the corresponding 
steady state particle distributions are presented in Figure 2. The spectra are primarily due 
to inverse Comptonization of soft photons modified by 7 — 7 absorption. Comparing the 
timescale for positron annihilation te+ ~ l/{nT(TTc) with that of inverse Compton cooling 
tic ~ 7/7/c reveals that 

!b^£5 3 X m^L,^^,R-'T;l (16) 
tic nT 

where I/^^ss, Ry and 1 are the soft photon luminosity, the size and soft photon tempera- 
ture normalized to 10'^^ ergs s"-*^, 10'^ cm and 1 keV respectively. 6 = max(l, -^^f^) where 
< Cs >~ 3A;T is the average soft photon energy. This implies that positron annihilation 
and the consequent annihilation emission is generally negligible compared to inverse Comp- 
tonization. In other words, positrons produced by p —p or 7 — 7 interactions cool by inverse 
Comptonization before annihilating. The luminosity due to p — p interactions is 

Lp_p ^ 8 X 10-^'^ergs s'^fr^Rr (17) 

The high energy photons produced by the inverse Compton process, pair produce by inter- 
acting with the soft photons. This results in a spectral cut off around ^ mgC^/S/cT, ~ 80 
McV. For high luminosities, the number density of 7-ray photons increases and pairs arc 
also produced by self interaction of the 7-ray photons. The importance of this process may 
be characterized by the compactness parameter. 
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such that when Ip^p 1, a luminosity dependent spectral break appears at an energy > 1 
MeV. The three spectra in Figure 1, correspond to increasing /p_p of 0.8, 4 and 8 and show 
spectral breaks at different energies. The spectra with these breaks is qualitatively similar 
to that obtained in earlier works (Lightman & Zdziarski 1987) where the soft radiation was 
assumed to be UV photons. 

To check a posteriori that other reactions apart from the p — p is not important in such 
scenarios, we compute the pair production rate due to the p— 7 reaction using the photon and 
non-thermal proton densities corresponding to / = 0.05 and other parameters given in Figure 
1. The reaction-rate expression for the processes was taken from (Chodorowski et al. 1992). 
The computed pair production rate ~ 1.7 x 10^^ cm~^/s is nearly two orders of magnitude 
smaller than that due to p — p reaction Ri 1.2 x 10^^ cm~^/s. However, this result depends on 
the proton energy index a and for flatter distributions, p — 7 reactions may be important. 

The resultant spectra generally have a broad feature around 1 — 50 MeV and hence 
may be detected by INTEGRAL and from future observations by GLAST. Presently, OSSE 
data when combined with simultaneous low energy data, can provide upper limits on the 
non-thermal fraction /. As an illustration, we consider the case of GRS1915-I-105 which was 
observed in a soft state simultaneously by OSSE and RXTE on April 21 1999 (OSSE VP 
813). (Zdziarski eA, al. 2001) fitted this data to the hybrid EQPAIR model and obtained the 
best fit black body temperature = 1.35 keV, the ion optical depth r = 4.4 and the luminosity 
of the soft photon source at a assumed distance of 12.5 kpc to be 1.3 x 10'^*^ ergs . On 
the assumption that these values may be approximated to be the soft photon temperature, 
optical depth and soft photon luminosity of this model, we compute the spectra for R = 10^ 
cm, a — 2.5 and for different values of /. Comparison of the observed OSSE flux at 500 keV 
{EFe ~ 0.1 keV cm~^ s~^) with the computed ones, provides an upper limit on the non- 
thermal fraction / < 0.05. The corresponding hmit on p — p luminosity is Lp_p < 8 x 10^^ 
ergs s~^. On the other hand, the power injected in non-thermal electrons as inferred from the 
EQPAIR model fit is 5 x 10^^ ergs s~^ (Zdziarski et al. 2001). Thus, on the basis of this 
simplistic and preliminary analysis, it would seem that the unknown accelerating mechanism 
preferentially energizes the electrons. However, it should be emphasized that an affirmative 
statement can only be made after a more careful analysis. 

As mentioned earlier, the positron annihilation time-scale is longer than the inverse 
Compton one and hence the spectra shown in Figure 1 do not show any detectable annihi- 
lation lines. For low energy positrons 7 fa 1 the ratio of the time-scales is 

^ 177L,,38i?7-' (19) 

tic 

which is still larger than unity. Hence the positrons may cool and thermalize before annihi- 
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lating and give rise to narrow {a ~ IkeV) annihilation line. The high resolution spectroscopy 
of INTEGRAL may be able to detect such a line, which would be a self-consistency check 
and may put constraints on the size of the system. However, gravitational red-shift and 
broadening due to Keplerian motion of the disk may inhibit such a detection. 

4. Summary and Conclusions 

It is shown that the presence of non-thermal protons in the inner regions of an accretion 
disk around a black hole, may have a detectable high energy spectral signature. The non- 
thermal protons produce electrons, positrons and 7-rays hj p — p collisions. The ambient 
copious soft photons are Compton up-scattered by these pairs to 7-rays. Subsequently, 
these photons produce more e^e+ pairs and an equilibrium is reached with the steady state 
emergent spectrum having a broad feature around 1 — 50 MeV. This spectral feature may be 
detected by future observations from INTEGRAL or GLAST. It is illustrated that presently 
available high energy observations of GRS 1915-1-105 by OSSE, impose an upper limit on 
the fraction of non-thermal protons to be < 0.05. 

The detection of such a. p — p high energy spectral signature will open a new window 
on our understanding of the acceleration mechanism and environment around black hole 
systems. It would also warrant a more careful and sophisticated treatment of the spectral 
computation than what has been done here. In particular, the specific geometry of the 
plasma (e.g. spherical cloud, corona, flare) has to be considered. The effect oip—p reactions 
between non-thermal protons and the cold disk ones need to be taken into account and a 
more realistic treatment of the soft photon source (which is probably the cold disk) needs to 
be undertaken. The escape of rclativistic protons (perhaps along open magnetic field lines) 
may also need to be considered. As shown in this work, the non-detection of this spectral 
feature will put constraints on the fraction of non-thermal protons existing in a source. This 
in turn will put strong constraints on the acceleration mechanism, since that would mean 
that the unknown mechanism preferentially accelerates electrons. 

The authors thank P. Subramanian for useful discussions. 
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Fig. 1. — Computed spectra for three different values of the non-thermal proton / = 
0.1,0.05,0.01 labeled as 1,2 and 3 respectively. The other parameters are R — 10^ cm, 
Ls = 10^^ ergs s-^ kT^ = 1 keV, r = 2 and a = 2.5 
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Fig. 2. — The steady state positron (bold) and electron (dashed) number densities in the 
plasma corresponding to the spectra in Figure 1 



